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Abstract 

We perform a 3D reduction of the two-fermion Bethe-Salpeter equation 
based on Sazdjian's explicitly covariant propagator, combined with a co- 
variant substitute of the projector on the positive-energy free states. We 
use this combination in the two fermions in an external potential and in the 
three-fermion problems. The covariance of the two-fermion propagators 
insures the covariance of the two-body equations obtained by switching 
off the external potential, or by switching off all interactions between any 
pair of two fermions and the third one, even if the series giving the 3D 
potential is limited to the Born term or more generally truncated. The co- 
variant substitute of the positive energy projector preserves the equations 
against continuum dissolution without breaking the covariance. 

PACS ll.10.Qr Rclativistic wave equations. 

PACS 11. 10. St Bound and unstable states; Bethe-Salpeter equations. 
PACS 12.20.Ds Specific calculations and limits of quantum electrodynamics. 
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1 Introduction 

The Bethe-Salpeter equation jl], ^ is the usual tool for the study of relativistic 
bound states. The principal difficulty in the treatment of this equation comes 
from the existence of unphysical relative time variables. The most usual practice 
in the two-fermion problem consists in eliminating the relative time variable 
(3D reduction). This 3D reduction is most often based on the replacement of 
the free Green function by an expression combining a delta fixing the relative 
energy and a 3D propagator. The exact equivalence (in what concerns the 
physically measurable quantities of the pure two-fermion problem) with the 
original Bethe-Salpeter equation can be obtained by recuperating the difference 
with the original free Green function in a series of correction terms to the 3D 
potential. 

The 3D reduction of the two-fermion Bcthe-Salpeter equation has been per- 
formed by many authors [3-19]. All methods are theoretically equivalent at 
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the limit of all correction terms included. Beyond the two-fermion problem, 
we have the cases of two fermions in an external potential ( 2^— body problem) 
and of three fermions (3-body problem). In these cases we meet new difficul- 
ties which do not appear or can be easily solved in the two-fermion case: the 
Lorentz invariance / cluster separability problem and the continuum dissolution 
problem. 

— Lorentz invariance and cluster separability. It is always possible to render 
an equation Lorentz invariant by working in the general rest frame (center of 
mass reference frame) and by building invariants with the total 4-momentum 
vector. The real difficulty appears when we start with a three-fermion equation 
and "switch off' all interactions between fcrmion 3 (for example) and the other 
ones. We should get a free equation for fermion 3 and an acceptable equation 
for the (12) system (cluster separability). In particular, this implies that the 
equation for the (12) cluster must not refer to the global center of mass frame 
anymore, as the momentum of fermion 3 enters in the definition of this frame. 
This equation must therefore be invariant, explicitly or implicitly (i.e. after 
rearrangements) . 

— Solution of the continuum dissolution problem. In the relativistic equa- 
tions for several relativistic particles, the physical bound states are degenerate 
with a continuum of states combining asymptotically free particles with oppo- 
site energy signs. This often neglected fact forbids the building of normalizable 
solutions in the N > 2-body problem (including the two-body plus potential 
problem). The usual solution consists in including positive-energy projection 
operators into the zero-order propagator [20-26]. The modified equations must 
of course continue to satisfy the Lorentz invariance / cluster separability re- 
quirement. 

In our preceding works on the two-particles in an external potential problem 
[ p3| , and, more recently, on the three-fermion problem [Efl, our approach was 
based on the cluster separability requirement. In the three-fermion problem 
the 3D potential was obtained by adding the three 3D potentials obtained by 
the 3D reduction of the three two-fermion equations in the three-fermion rest 
frame. Each two-fermion potential depending in general of the total energy of 
the subsystem, this two-fermion energy was taken as the three-fermion total 
energy, minus the free Dirac hamiltonian of the spectator fermion. The contin- 
uum dissolution was avoided in both cases by choosing a two-fermion reduction 
based on a 3D propagator containing a projector on the positive eigenvalues 
eigenstates of the Dirac's free hamiltonian. With these choices the three two- 
cluster limits are exact: switching off the two interaction terms with a spectator 
fermion gives the exact 3D equation which would be obtained by the reduction 
of the two-body Bethe-Salpeter equation. From this exact 3D equation for two 
fermions it is possible to go back to the original two-fermion Bethe-Salpeter 
equation and to perform the 3D reduction again in another reference frame, so 
that we can consider this two-fermion equation as implicitly covariant. This 
covariance is however quite implicit and global: neither the approached propa- 
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gator nor the positive-energy projector being covariant, the terms of the series 
giving the 3D propagator are not individually covariant, so that the implicit 
covariance becomes an approximated covariance when the series giving the 3D 
potential is truncated, for example by keeping only the Born term. 

In the present work we combine Sazdjian's explicitly covariant propagator, 
based on the second-order equations, with a covariant form of the positive- 
energy projectors. When used in the 3D reduction of a two-fcrmion Bcthc- 
Salpeter equation, this combination leads to a 3D potential given by a series in 
which each term is individually covariant. This potential can then be used in the 
two-fermion in an external potential problem and in the three-fermion problem, 
leading to continuum dissolution-free equations which can be truncated (for 
example by keeping only the Born terms) without losing the covariance of their 
cluster-separated limits. 

Although our 3D equations have exact two-cluster limits, they are themselves 
approximations. In the 2-i — body problem, one should take into account the 
modifications brought by the external potential to the fermion propagators also 
inside the Bethe-Salpeter irreducible kernel. In the three-body problem one 
has to take into account the irreducible three-body terms at the Bethe-Salpeter 
equation level and also the three-body terms generated at the 3D level by the 
reduction. We tried to do that recently j26j. 

In section 2 we build a covariant propagator by combining Sazdjian's covariant 
propagator with a covariant form of positive-energy projector. We use it to 
perform a 3D reduction of the two-fermion Bethe-Salpeter equation and study 
the resulting 3D equation (one-body limit, transition matrix elements, sym- 
metrization of the potential). We present also a "covariant Salpeter propagator" , 
slightly more complicated but preserving the particle-antiparticle symmetry. In 
the two next sections (3 and 4) we exploit these results by using our 3D po- 
tentials in the two fermions in an external potential problem (with examples) 
and in the three-fermion problem (with examples and with a calculation of the 
two-body limits). Section 5 is devoted to conclusions. 

2 Two fermions. 
2.1 Notations. 

We shall write the Bethe-Salpeter equation for the bound states of two fermions 
1 as 

$ = G°K§, (1) 

where $ is the Bethe-Salpeter amplitude, function of the positions x%,X2 or 
of the momenta p\,P2 of the fermions, according to the representation chosen. 
The operator K is the Bethe-Salpeter kernel, given by the sum of the irreducible 
two-fermion Feynman graphs. The operator G° is the free propagator, given by 
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the product G[G 2 of the two individual fermion propagators 

G° = ^--rA = 2 PM+ 2 k : . A (2) 

Pio - hi + iehi p\ - mf + it 

where the hi are the Dirac free hamiltonians 

hi = di .pi + Pi mi (i = 1, 2). (3) 

We shall define the total (or external, CM, global) and relative (or internal) 
variables: 

X = -(an + x 2 ) , P = Pi+P2, (4) 

X = Xx-X 2 , P=^(Pl-P2)- (5) 

and give a name to the corresponding combinations of the free hamiltonians: 
S = hi + h 2 , s = ^(hi-h 2 ). (6) 

We know that, at the no-interaction limit, we shall have to get a pair of free 
Dirac equations: 

(pio - h)^ = 0, (p 20 - h 2 )* = 0, (7) 
where VP depends on x%,X2. Let us also write their iterated versions 

(p\ - = 0, 04 - £f )M> = (8) 

with 

£ i = V ^f = (p? + m2)i (9) 

Interesting combinations can be obtained from the sum and differences of the 
equations (Q) or of the iterated equations (||): 

(P - S)V = 0, Cpq - s)* = 0, (10) 

H°* = 0, ( Po -n)* = (11) 

with 

= 2[{p\ - ml) + (pi - m 2 2 )] P0=fl = P 2 - 2{El + E 2 2 ) + V, (12) 

m = ^r(i?i 2 - £§) = ^(A? - hi) = f. (13) 
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2.2 Sazdjian's explicitly covariant propagator. 

The free propagator G° can be written 

G° = G 02 (p 10 + h 1 )( P20 + h 2 )(3 1 f] 2 , (14) 
where G 02 is the second-order (spinless) free propagator: 

G = (pi -E 2 +ie)(p 2 Q -E 2 +ie) ' (15) 

In a two-boson Bethe-Salpeter equation the free propagator G 02 could be ap- 
proached with the product G S2 of a constraint S(po — /i) fixing the relative 
energy and forcing a solution of the free difference equation (second equation 
(PH) by a 3D propagator which could be the integral of G 02 with respect to 
the relative energy po : 

G S2 = %0-M) Jd Po G 02 ( Po ) = -4i7r%o-M)^o ( 16 ) 

with 

<r=-(a 1+ a 2 ), n = -^~, ^ = ^-. (17) 

On the mass shell Po = hi+h 2 , the operator <7j is the sign hi/Ei of the energy 
of the fermion i. The presence of a annihilates thus the residues of the poles 
at P = ±(Ei—E 2 ), so that G S2 can be written 

G S2 = -2m 5(po-n) El±^l 1 . (i 8 ) 

U 2E X E 2 P 2 - (Pi + E 2 f + ie 1 1 

The propagator ( |l6| ) could easily be written in covariant form if it did not 
contain the operator a. If we replace a by e(Po) an d add the spinor part of 
(pT|), we obtain Sazdjian's propagator jl2| : 

G sz = 5( Po -[i) (pw + /n)(p 20 + (19) 

which can also be written in covariant form 

G^ = -M* J (P ■ p - ^j^i ) (Pl 2 • 71 + mi) ^ 2 • V m2) (20) 
v 2 7 p 2 +p 2 - (to 2 + to 2 ) + ie v ' 

or in a form which exhibits the pole in P — S : 

""■ij^EiiW*-"** <21> 
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2.3 Covariant positive-energy projectors. 

The combination of two fermions with opposite energy signs leads to continuum 
dissolution when an external potential or a third particle is added. The operator 
a would have kill the unwanted poles of Sazdjian's propagator, but we could 
not re-introduce it without losing the covariance. We have thus to search for an 
invariant substitute of a. A Lorentz invariant operator forcing a positive energy 
for fermion i could be 

Oi = e{p i0 )6{pD=6{p i0 -\p i \) (22) 

as the sign of the energy is invariant for positive squared 4-momenta. Taking 
the product of two such operators (one for each fermion) and replacing p by 
\x gives 

' p o P1-P2 i ~ , . rn\ - rr. 

T ~~ 2pJ ' P2 ' ~^2P^~) ' (23) 

Let us assume that the heaviest fermion is fermion 2. We can then write (|23| 

as 

P ° < .*|W$ + 4#-|AI) (24) 



2 2Pn V 2 IP, 



with 



Q2 = V P 2 2 + - mf . (25) 



The sum of the arguments of the two 9 is -Fo — I Pi I — I _P2 1 - The argument of the 
first 9 is (Pq — \pi\ + q2)(Po — \pi\ — <?2)/2-Po) which implies that Po must be 
outside the interval ([ p\ \ — q2, \ pi \ +92)- The combination of these two results 
implies Po > |pi|+<?2- The second brings no supplementary restriction, its 
argument being always positive when Pq > |pi|+<72- We can thus write 



9 = 6(P -\pi\-^p$+m 2 2 -ml). (26) 

We see that our projection operator introduces a cutoff on the high \pi\ for a 
given P . 

It is interesting to write 9 also in terms of the hi , for comparison with the 
A ++ =9{hi)9{h2) projector widely used by us |2^, ^6| and by others, and also 
for future use in the two-fermion plus potential problem. We get 

9i = 6( Pl0 )9(p 2 t ) = 9( Pl0 ) 9{p% -hj + m}) (27) 

so that 

9 = [9(p w ) 9(p 2 w -h\+m\) 9(p 2Q ) 9(p 2 2Q - h 2 + m 2 )] . (28) 
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The po = fi constraint gives 



P h 2 -hi P h 2 -h 2 

^ = f + ^T< ^-f + ^T (29) 

and also p 2 Q — h 2 = p\ — h\ so that the last 9 which concerns the heaviest 
fermion can be omitted. For a given bound state energy Po, 9 eliminates the 
mixed energy-sign continuum which would share this energy. When Po = S, 
( p8| ) reduces indeed to 9(h\)9(h2)- In fig. 1, we draw a map of the 9 = 1 region 
in the (hi, /12) plane, for Po fixed to 2 (in arbitrary units) and different values 
of the lowest mass mi : 0, 0.5, 1 (solid line), 1.5, 2, 2.5 and 3 (while m-z can 
take any value larger or equal to mi). We draw only the first quadrant, the 
other ones being the mirror images of the first one. The 9 = 1 region is bounded 
by the two Pio = curves, independent of the masses, and by the p 2 = curve 
corresponding to the lowest mass mi. In each quadrant, this region is divided 
into four parts, corresponding to the signs of h 2 —m 2 and — m| (as an 
exemple, we indicate this partition in the mi = 1, 7712 = 1.1 case, for which Po 
is thus 0.1 below the threshold 2.1). When no external potential is present, both 
expressions are positive definite, the negative values being for the bound states 
of the corresponding fermion in an external potential. 

While the widely used A ++ = 0(/i 1 )0(/i 2 ) noncovariant projector simply se- 
lects the first quadrant, the 9 = 1 region is divided into four symmetric parts, 
one for each quadrant, but there is a cutoff on the high-|/ij| values. For moder- 
ately relativistic systems, the important regions are near the (mi, ^12) point. 

We shall thus finally approach the free propagator by the product of the 
covariant Sazdjian propagator with our covariant projector: 

G s = § G sz ( 30 ) 

2.4 3D reduction of the two-fermion Bethe-Salpeter equa- 
tion. 

We shall write the free propagator as the sum of the zero-order propagator, plus 
a remainder: 

G° = G s + G R . (31) 
The Bethe-Salpeter equation becomes then the inhomogeneous equation 

$ = G°K& = (G 5 + G R )K$ = * + G R K (32) 

with 

* = G S K$ (= G 5 (G r 1$ )- (33) 

Solving (formally) the inhomogeneous equation ( [32] ) with respect to $ and 
putting the result into (J33]), we get 

* = G 5 K(1 - G R K)- 1 ^ = G S K T ^> (34) 
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where 

K T = K(l - G^K)- 1 = K + KG R K + ... = (1 - KG R )~ 1 K (35) 

obeys 

K T = K + KG R K T = K + K T G R K. (36) 

The reduction series ([35]) re-introduces in fact the reducible graphs into the 
Bethe-Salpeter kernel, but with G° replaced by G R . Equation (^) is a 3D 
equivalent of the Bethe-Salpeter equation. It depends on the choice of G 5 . 
The relative energy dependence of eq. ([mJ) can be easily eliminated: 

y = S( P o-v)ip (37) 

and if) satisfies the equation 
with 

/W^/i) = f dp' Q dp S(p' -^p 1 p 2 K T (p' QlPo )S(p -^. (39) 



Note that we write (j?o,Po) but {fi,fx), as we keep fi in operator form. This 
operator is diagonal in the spatial momentum space. The eigenvalue will depend 
on the position of fi in the formula: the eigenvalue of the first /i in (|3^) will 
be built with the final momenta and that of the last /i will be built with the 
initial momenta. Defining 

V = -2in 9 (3i(3 2 K T {n, fj,) 9 (40) 

we get 

^= - --±--Vi>. (41) 

V Po-S + ieP 2P v K ' 

We included a final 9 in the definition of V, as the wave function at right is 
already projected, and we replaced |Po| by Pq, which is always positive in the 
9 = 1 region. 

The inversion of the reduction is given by 

$ = {1-G R K)- 1 f = (l + G R K T )^ = (1+G°K T -G S K T )^ = G°K T ^ 

(42) 

or, explicitating the relative energy variable 

$(p> ) = G Q (p' )K T (p' ^)^. (43) 
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2.5 Covariant Salpeter equation. 

A Bethe-Salpeter equation with an instantaneous kernel (i.e. a kernel indepen- 
dent of the relative energy) leads directly to Salpeter's 3D equation g] : 

*= P -S + ieP VS * (44) 

with 

V s = ^ilTTPifcKT 2 (45) 

r = 6(hi) 6{h 2 ) - 9{-h 1 ) 6(-h 2 ). (46) 

When the kernel is not instantaneous, we can make an expansion based on the 
approached propagator 

G 5S = 5(p Q -n) f d Po G°( Po ) = -2 t nS(p -^- T . Pifo (47) 

J P a - S + ieP 

in the same way as above. This operator r, which projects on the r 2 = 1 
subspace with a change of sign for the negative-energy part, is also a continuum 
dissolution-preserving operator. 

A "covariant Salpeter equation" can easily be obtained by combining Sazd- 
jian's covariant operator with a covariant substitue of r instead of a covariant 
substitute of 9(hx)9(h 2 ). We could indeed use 



(48) 

Po=P 



[0(pio) *(Pao) - 0(-pio) 6(-p2o)] HpD 0(pI 

Our 3D equation would then again be ([4l|), but with the potential 

V = —2in 9 s (3i(3 2 K T (fi, fi) (9 s ) 2 (49) 

where K T must now be built using 9 s instead of 9. 

As the second term of ([l6]) does not contribute much in practical calculations, 
it is often omitted. For the same reason, we shall omit the second term of (Esj) 
and use simply 9 in the present work. From a more fundamental point of view, 
however, the use of r or of 9 s preserves a particle-antiparticle symmetry which 
is an important feature of relativistic theories. 

2.6 One-body limit. 

If we go to the center of mass reference frame, write Po = rn 2 + W\ (for Po > ) 
and make m 2 — ► oo in the 3D equation (|4ll), we get the equation 

^ = Tf7 \— V oo V>, Voo = lim V. (50) 
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In |L7|], we showed that a 3D reduction performed with a propagator AG SZ , A 
being a projector, leads to a 3D potential given by 



V = A 



Po-S 



A. (51) 



We also showed that the limit of V (let us call it V c ) is given by the 
limit of the Born term (a Coulomb potential in QED). Since the limit of 9 is 
9(W\ — |f?i|), we have thus 

K» = $(W! - V c (l - V c ^j 6(Wi - \pi\). (52) 

Equation ( |50| ) is related to the equation 

^ = = \ - V C ^ C (53) 

which describes the light fermion 1 in the potential generated by the heavy 
fcrmion 2. We have 

iP = 0(Wi - \pi\)^ c . (54) 

Equation ( |5C| ) is thus the exact equation for the projection (|54|) provided the 
expansion of ( |52] ) is not truncated. If, however, the two-body potential is trun- 
cated to the Born term, its one-body limit will be the projection of V c instead 
of Voo- The introduction of the projector 9 has destroyed the exact one-body 
limit of the Born approximation, as it would also be the case with other contin- 
uum dissolution-preserving operators like A ++ . 

2.7 Transition matrix elements. 

The 3D off mass shell transition matrix element is 

T 3D = V + V V +... 

Po-S + iePo 2P 

= -2iir 9 p!(3 2 [K T + K T G S K T + ■ ■ ■] {ft, p) 9 (55) 
K T + K T G S K T + ■■■ = K {\-G R K)-' s -{\ - G s K(1~G R K)~ 1 )^ 1 

K(1-G R K-G 5 K)- 1 = K{1-G°K)- 1 = T (56) 

so that 

T 3D = -2iir 9 13x02 T{p, n) 9. (57) 

On the positive-energy mass shell Pq = E\ + 9 = 1 and T(fi,fi) is the 
physical on mass shell scattering matrix element of field theory. Although the 
Bethe-Salpeter equation (^) was written only for bound states, our final 3D 
equation can still be used for describing the scattering. 
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2.8 Symmetric form. 

In our 3D equation (^), written in the form 

(P -S)^ = ^-^-V^ (58) 

the interaction term (the operator before ip at right) is not symmetric and 
depends on Pq via the operator before V and also via V itself. This does not 
excludes a real energy spectrum, but this can not be seen directly on J58]). In 
the two-fermion framework, this is not a problem as there exist many ways of 
transforming ([58]) into a symmetric equation. We shall however use these two- 
fermion interaction terms in a three-fermion equation below (in section 4) and, 
in this new framework, their symmetrization will not be possible anymore (it 
would demand a different transformation for each two-fermion pair). In order 
to be sure to get a real energy spectrum from our three-fermion equation, we 
should thus symmetrize the interaction terms at the two-fermion level. This can 
be done by simply multiplying ( p8| ) by 1 + V/2P , which gives 

(P -S)iP = Uip, (59) 

In a 3D reduction based on the noncovariant projector (p^), we would have 
simply V s . The first-order energy shift — < V 2 /2Pq > corresponding to the 
difference U — V would then be included in the contribution of the first ladder 
term KG RS K, which enters in the definition of V s . 



3 Two fermions in an external potential. 
3.1 3D equation. 

The two-fermion in an external potential problem is already much more compli- 
cated than the pure two-fermion problem, although it exhibits some simplifying 
features when compared with the three-fermion problem. The principal new 
difficulty is the non-conservation of the total spatial momentum P. In the pure 
two-fermion case, this conservation law forbids the mixing of the physical bound 
states with the mixed-energy states (continuum dissolution) . When an external 
potential is present, it becomes possible to get any given energy in an infinity of 
ways by combining fermions with opposite free energy signs. The positive-energy 
bound states will then not be normalizable, being mixed with a continuum. 

The two-fermion plus external potential equations can easily be obtained 
with the simple generalizations |23| 

hi — > h ie = di .pi + PiTTii + Vi(xi,ji) (61) 
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where Vi is the external potential acting on fermion i. These substitutions must 
be done in the free Bethe-Salpeter propagator G° and in the manipulations 
leading to the 3D equation (this is an approximation, as we should also modify 
the fermion propagators inside the irreducible graphs, but this gives the exact 
cluster-separated limits). We get 



Po + S e 
2 Pa 



V, 



1>, 



(62) 



with 



S e = hi e + h 2e = hi + h% + Vi(xu"ji) + ^(£2,72), (63) 
V e = -2in§ e / dp dp S(p' Q - n^PxfoK]; (p' ,p )5(p - fj, e ) § e , (64) 



Me 



hi) 



|(^ + Me )0(^- Me 
K T = K 



'<[? 

KG^K - 



hie +"»?)■ 



Gf(Po) 



1 



1 



\Pv +Po- h le + ieh le ±P ~ Pa ~ h 2e + ieh 2e 
2iir P + S e 



P -S e + ieP 2P 



S{po~^e) I 



0103 



(65) 

(66) 
(67) 



(68) 



Besides the free-free continuum, the (/ii e ,/i2e) spectrum (for which we can 
use figure 1) contains now bound-free or free-bound combinations (lines) and 
bound-bound combinations (points). The equations are written in the labo- 
ratory reference frame, defined as the reference frame in which the external 
potential is static. When the mutual interaction is switched off, we get a pair 
of independent Dirac equations in the external potential. When the external 
potential is switched off, we get the equation of a system of two mutually in- 
teracting fcrmions. This last equation can be written in covariant form, the 
contribution of each term of (67) being separately invariant. We can keep only 
the Born term (only K in i^Jand only the ladder term in K ) without losing 
this covariance. 



3.2 Examples. 

3.2.1 Two electrons in the field of a nucleus. 

In this case, the external potentials are given by (in configuration space): 

-Ze 2 -Ze 2 



Vi 



V 2 = 



xi 



*2 



(69) 
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For the mutual interaction kernel we have to choose a gauge. In Feynman's 
gauge, the Born term is 

Ktf,p)= -1^^(71-72), k= P = P >. (70) 

In Coulomb's gauge, it is 

is! > \ 2ie2 r 1 f- - lrkj2-k 

K{P > P) = - M 1 ^"^ 171 ' 72 "^ 1 ' (71) 

This last gauge being not covariant, the space-time separation refers to the 
laboratory frame. This gauge is not a priori a good choice for our covariant 
calculation. In the two-fermion problem, however, it is well known that the 
best choice for the calculation of the energy spectrum is Coulomb's gauge in the 
center of mass reference frame. We could take this fact into account by writing 
( fri]) for a pure two-fermion system in the two-fermion center of mass reference 
frame P = 0, making everything covariant by using P : 

P\P2 -> ^ ' a-b^ ^ (°' & )' ( 72 ) 

and importing the result into the two-fermion in an external potential problem, 
where P is no more conserved. We could call this the "covariant Coulomb 
gauge" . 

We have now to compute V e by using (|^) , with = K given by ( |70| ) . 
The principal difficulty comes from the fact that the modified Dirac operators 
hie are diagonal neither in momentum space, nor in configuration space. If we 
expand the potential V e into the eigenfunctions of the hi e , it will be given by 
(in Feynman's gauge): 

V e (pl,Pi) = MiP^'M) < < I V e \"i > <(£i)V>i(P2) (73) 

with 

< u[ | V e | wt > = - ^ 6M) J dZp'^p'^pidZpi < (p[)^ 2 [p' 2 ) 

i — di '6l<i 



UJ i UJ-i 



, ,\ > , =^2 VV;iG»l)^ 2 (P2)0e(Wi) (74) 

[MeK) - Mujt)} 2 - {p -pY 

where ip Ui (J>i ) is an eigenstate of hi e with eigenvalue u>i . The sum in ([73]) 
includes also an integral on the continuum. This expression is well adapted to 
a perturbation calculation starting with an eigenstate of S e . The first-order 
energy shift will indeed be < u>i \ V e \ u>i >, in which 9 e (uJi) will be equal to 1, 
as Pq=oji+lj 2 . 
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3.2.2 Positronium in an external potential. 



For the positronium in an external potential, the mutual interaction kernel 
and potential in the Born approximation (neglecting the annihilation graph) 
are the same as for the two-electron pair above, with a change of sign. The 
external potentials could again be the potentials generated by a nucleus, and 
we would write again (69), with a change of sign for the positron. We could 
also consider non-central potentials (as V\ = &X\ 3 , V-z = — ex 23 for a constant 
electric field along the third axis). The representation of V e in the space of the 
eigenfunctions of the modified Dirac operators gives a meaning to our rather 
formal expression (|64|) , but it is not well adapted to this new framework, as the 
principal interaction is now the binding mutual interaction. We shall therefore 
isolate a zero-order Coulomb contribution from V e ■ 

V e (p-,fi ) = V C {p',p)+Y J i>»'M)^'M) < J i I (Ve~V C ) | Wi > TpXiPMM) 

(75) 

with 

<u' i \(v e -v c )\u i >= ^oM) I "(PpWicPpi^PitiMwM) 



1 — ai-d<2 



4 Three fermions. 



ipwtipi)^^) Oe(t^i) (76) 



4.1 3D equation. 

It is easy to build a three-fermion 3D equation with three two-fermion potentials 
Uij(Pijo) (we explicitate for a while the dependence in the total energy of the 
two fermions): 

(P -S)i> = [U 12 (Po-h 3 ) + U 23 (Po-hi) + U 31 (P ~h 2 )}^ (77) 

with S=hi + hi + h,3. We have replaced the arguments Pijo, by the operators 
Po—hk in order to get the exact cluster-separated limits. At the V23 = V31 = 
limit, for example, we get indeed 



P — hi — h2 — h 



Writing 

ip = 1P12 ipa, 

we get two independent equations: 
■Pi 20 —h\—h2 



1> = U 12 (Po-h 3 )iJj. 

Pf) = Pi 20 + P30 
"012 = C^12(Pl20) ^12, P3O03 = h 3 1p 3 



(78) 
(79) 

(80) 
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as Po — hs can be replaced by P0—P30 = Pi20- Our 3D equation (77) satisfies 
thus clearly the cluster separability requirement. Furthermore, the three cluster- 
separated limits are exactly the 3D equations we would get by 3D-reducing the 
corresponding two-fermion Bethe-Salpeter equations. This cluster separability 
is a property of the equation, or, if we prefer, of the full Green function. For a 
given scattering solution it is also possible to take the cluster-separated limit at 
fixed P. This is not possible for the bound state solutions. 

Our use of a covariant approached propagator results in an explicit covari- 
ance of the cluster-separated limit equations like (^0|): this equation can directly 
be written in covariant form and each term of the series giving the 3D potential 
is separately covariant. We could thus neglect all terms of K T but the Born 
term without losing the covariance of the cluster-separated limits. 

Let us recapitulate the elements of our 3D equation (Q) (we give only U\ 2 ) : 

U 12 V 12 {P -S, } + 4(Po _ hz? (81) 

V 12 = -2i7r0i2/?i/kJif&(/ii 2 ,/ii2)0i2 (82) 



U 12 



= 9(P -h 3 - I - yjft + m§ - m?) (83) 

^ = 2(ph^ {El - El) (84) 

= K 12 + K 12 G? 2 K 12 + ■ ■ • (85) 
1 



(pio -hi + iehi) (p 20 - h 2 + ieh 2 ) 

2Z7T P Q + S - 2/l 3 



Pa - S + i€ 2{P - h 3 ) ) 

We could keep only the Born terms in the Kj- without losing the covariance of 
the cluster-separated limit equations. 

4.2 Two-body limits. 

It is interesting to compare the two-body limits (say m 3 — » 00 ) of our three- 
fermion equations with our two-fermion in an external potential equations. 
Writing P a = m 3 + P 12 o and making m 3 — > 00, we get 

Pim1> = {hi + h 2 + Vx + V 2 + U l2 )^ (87) 

where C/ 12 is still given by (jsil), with P — h 3 replaced by Pi 2 q, while 



>(-Pi20 + h 2 + \pi\) wC 



V l = 6{P 120 -h 2 -\p x \)Vt 1- v p \ Iy Vt e(P 120 -h 2 -\pi\ 

Piw - hi - hi J 

(88) 



-1 
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v 2 = e(p lw - hl -\p 2 \)v 2 c (i- '(-^M + fr + lftl) ^! ecj^-^-iAi 

V ^120 — n 1 ~ n 2 J 

(89) 

We note the following differences with the two-fermion in an external potential 
equation: 

- In the mutual interaction term (which takes its symmetrized form), the 
free Dirac operators hi are not replaced by the operators h ie . 

- The external potentials are projected, as in the one-body limit (|52|). The 
energy argument W\ is here P\ 2 Q — h 2 (and vice- versa). If is energy- 
dependent (to consider more general cases than a pure Coulomb potential) its 
energy argument will also be replaced by P\ 2 o~h 2 . 

The difference between both approaches consists clearly in three-body terms. 
The cluster-separated limits are still exact, as they were before taking the m 3 — > 
oo limit. These cluster-separated limits appear when two of the three two- 
body interactions are " switched off' . When only the mutual interaction V\ 2 is 
"switched off' , we do not completely get two independent equations as in section 
3. Writing P 120 =pio+P20, we see that V\ will depend on p 2 o~h 2 and V 2 on 
Pio~ h\. It must however be noted that the external potentials obtained by the 
m 3 — > 00 limit are less general than the external potential we could consider in 
the formalism of section 3: here we expect central potentials decreasing with the 
distance. With such potentials, it is not possible to let fermion 1 (for example) 
go to infinity without having both V\ 2 and V31 going to zero. 

4.3 Examples. 

4.3.1 Two electrons and a nucleus. 

For two electrons and a nucleus of spin i, we have the Born terms 

— 2zc 1 

2i Ze 2 1 

K 23 (p' 23 ,p 2 3) = ro )2 , . (72-73) (91) 



-1 



^'w^^' (92) 

The 3D potentials are, in the momentum representation (we omit the momen- 
tum arguments of the Vij here and in the next subsection) : 

V 12 = ^9(P Q - h 3 -\p[\- \p' 2 \) [~ dl f 2 2 9 (P -h 3 - \pt\ - \p 2 \ ) (93) 

77 K-120 ^12 



-2ie 2 


1 




(2tt) 3 


kf 2 + 


ie 


2i Ze 2 


1 




(2tt) 3 


^23 "+ 


■ ie 


2iZe 2 


1 




(2tt) 3 


k 2 A- 


■ ie 



Ze 



2 



V23 ^2^ e (P ° ~ hl ~ '^ 2 ' ~ V +m 2 -m 2 ) 
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Q(P n - hi - Inol - i /w 2 + m 2 - rr 2 

^230 ^23 



(P - /ii - \p 2 \ - J Pi + m\ - ml ) (94) 



where fc 2 20 and fc| 30 must be taken at 

1951 

while V31 = V23 (1 <-> 2). We could also use the "invariant Coulomb gauge" 
(in fact, one gauge for each two-fermion subsystem). If we had started with 
a true three-fermion Bethe-Salpeter equation we could not use three different 
gauges simultaneously, but here our aim is only to get continuum dissolution-free 
equations with the correct cluster-separated limits. 

4.3.2 Three protons exchanging mesons. 

For three protons exchanging a vector, a scalar and a pseudoscalar meson, the 
contribution of the Born terms to V12 is 

Vi2 = ^0(P -h 3 -\ Pl \-\p 2 \) (l-a r a 2 ^ ) 

I K 120 K 12 m p V 



glPlfc ffj/3l/327l5725 



fl(fl,-ft 3 -|ft|-|ft|), (96) 



v 120 ^12 "V ^120 "-12 "V 
,2 



where where fc 2 20 must again be replaced by (|95|), and similarly for V23 and 

I 31. 



5 Conclusions. 

Using the explicitly covariant Sazdjian propagator combined with our covari- 
ant substitute of the positive-energy projector, we built a 3D equation for two 
fermions in an external potential and for three fermions. The two-fermion equa- 
tions obtained at the cluster separated limits are these which would be obtained 
by a 3D reduction of the two-fermion Bethe-Salpeter equation. In our preceding 
works this implied that these cluster-separated limits were implicitly covariant 
if not truncated, even if they remained written in a specific reference frame 
(the laboratory frame or the three-fermion rest frame). In the present work, 
however, they are explicitly covariant, which implies that the covariance will 
survive the possible truncations of the 3D potentials. The invariant substitute 
of the positive energy projector we use brings a cutoff on the spatial momenta, 
which could be welcome when insuring the existence of some integrals. This 
cutoff appear for rather low, but nevertheless truly relativistic values. It is not 
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by itself an approximation: it reflects the choice to work with a given projection 
of the two-fermion equation. 

The 3D equations we write in this work are more complicated that the 
corresponding equations in our preceding works. Keeping in mind that in both 
cases the 3D potentials are given by series to be truncated anyway, we hope 
that the Lorentz invariance / cluster separability properties of the truncations 
of the present equations will lead to a better approach of the real physics in the 
first terms. 
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Figure 1. The 9=1 region. 




